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We report a search for the rare decays D^K'^ and B+ D^K*^ in an event sample of 

approximately 465 million BB pairs collected with the BaBar detector at the PEP-II asymmetric- 
energy e+e^ collider at SLAC. We find no significant evidence for either mode and we set 
90% probability upper Hmits on the branching fractions of BF{B+ D^K'^) < 2.9 x 10"^ and 
BF{B+ -> D+K*^) < 3.0 X lO"*^ [|l]]. 
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Figure 1: Annihilation diagram for the decay B^^D^K^*^^ (left) and hadron-level diagram (right) for the 
rescattering contribution to B+ — >^ D^K^*^'^ via B+ — )• DsTt'^. 

1. Introduction 

Charged B meson decays like —5- D^K^*^^ are dominated by weak annihilation diagrams, 
for which no rehable estimates for the decay rates exist because of soft gluons exchange. In particu- 
lar annihilation amplitudes cannot be evaluated with the commonly-used factorization approach [Q]. 
Such annihilation amplitudes are suppressed by where A is the sine of the Cabibbo angle |^ |3|]. 
So far, no pure annihilation hadronic diagram has been observed, and such amplitudes are usually 
neglected in the measurement of Vut- Their branching fractions could be enhanced by so-called 
rescattering effects (see Fig. [T]), up to [||], rendering the rate comparable to the isospin-related 
B° — ;> D^K^*^^ decay rate of approximately 5 x 10"''. 

None of the modes studied in this note has been observed so far, and a 90% confidence level up- 
per limit on the branching fraction )-D+/r°) < 5 x 10"^ has been established by BaBar 
No study of B+ — )• D^K*^ has previously been published. The results presented here are obtained 
with 426 fb"' of data collected at the T(45) resonance with the BaBar detector at the PEP-II asym- 
metric e+e" collider corresponding to 465 x 10^ BB pairs (Ngg). An additional 44.4 fb"' of 
data ("off-resonance") collected at a center-of-mass (CM) energy 40 MeV below the T{4S) reso- 
nance is used to study backgrounds from e^e^ qq {q = u, d, s, or c) processes, which we refer 
to as continuum events. The BaBar detector is described in detail elsewhere 

2. Event Reconstruction and Selection 

The D+ mesons are reconstructed in the modes — K~n^7l^ (Ktitt), D+ — Ks7r~^ (KsTt), 
D+ K-n+n+n° (Knnn'^) and D+ Ksn+n'^ {KsTfnP) for the decay channels B+ D+K^ 
(DK). Only the first two modes are used for the B+ — > D^K*'^ decay channel (DK*). The event 
selections are optimized by maximizing S/\/S + B, where S and B are the expected signal and back- 
ground yields, using Monte Carlo (MC) simulations and off-resonance data. The signal branching 
fraction is taken to be 5 x 10"^. 

The charged kaons are required to satisfy kaon identification criteria obtained from the combi- 
nation of information from the Cherenkov light and the tracking detectors. Kaons and pions must 
satisfy pK > 200 MeV/c and pjc> 150 MeV/c, where p is the momentum in lab frame. The invari- 
ant mass of the candidates is required to stand within 10 to 22 MeV/c^ (depending on the chan- 
nel) of the nominal mass [Q]. The Ks candidates are reconstructed from K^n^ pairs with invariant 
mass within 5 to 7 MeV/c^ of the nominal ^5 mass [|8]]. We define a^^.(B+) as the angle between 
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the momentum vector of the Ks candidate and the vector connecting the and Ks decay vertices. 
The prompt Ks candidates from the — )• D^Ks decay must fulfill ln(l — cos aKs{B^)) < — 8 and 
ln(l — cosajfj(D+)) < —6, where aKs{D^) is defined in a similar way. The tt" candidates are re- 
constructed from photon pairs 77 with invariant mass m(77) within 10 to 12 MeV/c^ of the nominal 
TT*' mass [|]. These pairs must satisfy £"(7) > 70 MeV, £"(77) > 200 MeV, PcMiYY) > 400 MeV, 
where E and Pcm are respectively the energy and the momentum in the CM frame. The K*^ can- 
didates are reconstructed in K*^ — K^n^ with the invariant mass liying within 40 MeV/c^ of the 
nominal K*^ mass |j8|]. We define % as the angle between the direction of flight of the charged 
K and the direction of flight of the B in the K*^ rest frame, and require | cos 0h| > 0.5. The 
candidates are reconstructed by combining one D+ and one Ks or K*^ candidate, constraining them 
to originate from a common vertex. We define Ob as the B polar angle with respect to the beam 
axis in the CM frame, and require | cos Ob] to be smaller than 0.76 to 0.86 depending on the chan- 
nels. Using the precise knowledge of the e^e^ beams energies and the energy conservation in the 
two-body decay T{4S) — )• BB, we define the beam-energy substituted mass niEs and the energy 
difference AE: 



where E and p are energy and momentum. We retain candidates with \AE\ value smaller than 
19 to 25 MeV and tues in the range [5.20,5.29] GeV/c^. Multiple B candidates are eliminated 
with selections on D+ mass or AE" distribution. The dominant background comes from continuum 
events, characterized by a jet-like topology, which can be described with these variables defined in 
the CM: the cosine of the angle between the B thrust axis and the thrust axis of all the other tracks 
and energy deposits of the event, where the thrust axis is defined as the direction that maximizes 
the sum of the longitudinal momenta of all the particles, the event shape moments Lq = ^, pi, and 
L2 = ^;;?,| cos where the index / runs over all tracks and energy deposits in the rest of the 
event; p, is the momentum and 0, is the angle to the B thrust axis. We also use \At\, the absolute 
value of the time interval between the two B decays ||5|]. These four variables are combined in a 
Fisher discriminant F [||], whose coefficients are determined with samples of simulated signal and 
continuum events, and validated using off-resonance data. For the Knn mode, events are classified 
according to their flavor-tagging category (lepton, kaon or other) and fitted simultaneously. The 
BB background is divided into two components according to their distribution in the signal region: 
non-peaking and peaking. The peaking backgrounds are rejected using the Ks helicity angle Ok^ 
with \cos{dKs)\ > 0.8 or 0.9 depending on the channel. Based on MC studies, atmost one BB 
peaking background event per mode is expected in the signal region. The charmless background is 
evaluated from data using the D+ sidebands and found to be negligible. 

3. Fit Procedure 

The signal and background yields are extracted with an unbinned maximum likelihood fit of 
niEs and F, assuming from simulation studies the correlations between mes and F to be negligible. 
For niEs the signal is modeled with a Gaussian function, the continuum and non-peaking BB back- 
ground are described by two ARGUS functions JlO|]: A{x) = — {x/xqY •exp(c(l — (x/xq)^)), 
where xq is the maximum value of x and c accounts for the shape of the distribution and are de- 
termined from data for the continuum. All other PDF parameters are derived from the simulated 
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events. The peaking BB background is modeled with a Crystal Ball function [|n|] which is a Gaus- 
sian modified to include a power-law tail. The peaking background yield is fixed from the PDG 
branching fractions [||]. The signal yield determined by the fit (Ngig) is used to calculate the branch- 
ing fraction (BF): BF = N^ig/ {Nb+ • Ssig • BF^et), where Nb+ is the total number of charged B mesons 
in the data sample, BFgec is the BF is of the secondary decay channels of the D and Ks, and Esig 
is the signal reconstruction efficiency measured in MC. The fit procedure is vaUdated using toy 
MC studies and no biases of the fit model were found. The fit model was tested using full MC 
sample with and without signal events. The results of the fit to the data are reported in Table ^ 
for each D channel. The background yields are close to the expectations and the errors obtained 
on the branching fractions are in good agreement with the values found with the toy study. The 
leading contribution is obtained from the Knn mode. The Fig. ^ gives the fit projection for thes, 
after requiring F > 0, to visually enhance any possible signal. 

Table 1: Branching fraction (BF) measured in units of 10^^ with statistical und systematic uncertainties for 
each channel. A^, are the yields of the fitted species. 



Decay mode 








BF 


B+ D+K^ 


Knn 




70 ±27 


2690 ±57 


-4.2+2;4(stat.)+[;'(syst.) 






111±51 


65 16 ±94 


20+2^'(stat.)+[[:3(syst.) 


KsTt 




20 ±14 


381±23 


0.7^[^(stat.)+^:2(syst.) 






36 ±22 


1270±41 


-14^^:2(stat.)+?fg(syst.) 


B+ D+K*^^ 


KitK 


-15.6+^:1 


463 ±63 


6338 ±98 


-5.0^|?(stat.)+|i(syst.) 


KsK 




35 ±15 


547 ±27 


-33^10i2(stat.)+«;4(syst.) 



4. Systematic uncertainties 

The uncertainties on the PDF parameterizations is evaluated by repeating the fit varying the 
MC-obtained PDF parameters within their statistical errors, taking into account correlations among 
the parameters. Differences between the data and MC for the signal PDF shapes are investigated 
using data control samples B^ — D^n and B^ — D^p^. The uncertainty on the continuum back- 
ground shape is estimated using off -resonance data instead of continuum MC. The uncertainty on 
the PDF of the non-peaking BB background is measured by leaving its parameters free in the fit 
and taking the difference from the nominal fits as uncertainty. We also considered the uncertainty 
on signal efficiency due to limited MC statistics. Uncertainties on MC-data differences in tracking 
efficiency, Ks and 71° reconstruction and charged-kaon identification, are estimated by comparing 
data and simulation performance in control samples. The uncertainty on peaking background are 
estimated by repeating the fit varying the event yields within their statistical errors. The uncertain- 
ties on the branching fractions of the sub-decay modes are also taken into account. The uncertainty 
on Ngg has a negligible effect on the total error. The uncertainties are included by convolving the 
individual fit likelihoods with Gaussians of width equal to the systematic uncertainty. The total 
systematic uncertainties on the BF are given in the Table |l| for each channel. 
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Figure 2: From top left to bottom right: niEs projection with F > for Knn, Krnzn^, K^K, and Ksnn'^ 
for B+ D+Ks and /:;r;r and Ksn for B+ ^ D+K*^. Data are black dots with error bars, the different 
fit components are : signal (black curve), non-peaking BB (green), continuum (magenta) and BB peaking 
background (red) and the total pdf (blue). 



5. Results for Branching Fractions 

The individual likelihoods for each mode are finally combined to give the average BF's, which 
are compatible with zero. We then quote an upper limits at 90% probability using a Bayesian 
approach with a flat prior for the BF: 

BF{B+ D+K^) < 2.9 x 10"^ BF{B+ D+K*^) < 3.0 x 10^^ 
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